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We consider a theoretical problem of the upper critical magnetic field parallel to a conducting 
axis of a quasi-one-dimensional layered superconductor. We show that the orbital effects against 
superconductivity in a magnetic field are capable of destroying the superconducting phase at low 
temperatures if the interplane distance is less than the corresponding coherence length. Applica- 
tions of our results to the recent experiments, performed in the superconductor Lio.gMoeOiy [J.-F. 
Mercure et al., Phys. Rev. Lett. 108, 187003 (2012)], provide strong arguments in favor of a triplet 
superconducting pairing in this quasi-one-dimensional layered conductor. 
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Quasi-one-dimensional (Q1D) superconductors have 
been intensively studied since the discovery of supercon- 
ductivity in the organic superconductor (TMTSF)2PF6 
[1,2]. Early experiments [2-5], performed on the Q1D 
superconductors (TMTSF) 2 X (X=PF 6 and C10 4 ), pro- 
vided hints of their unconventional nature. In addi- 
tion to possible triplet pairing [6-10], it was proposed 
that Q1D superconductors can demonstrate such un- 
usual phenomena as the reentrant superconductivity 
[7-10], Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase 
[2,7-15], and hidden reentrant superconductivity [15]. 
Although triplet superconducting pairing was consid- 
ered as the most probable mechanism for many years, 
more recent experiments [12,16] and theories [15,17] are 
in favor of a d-wave like superconducting pairing in 
the (TMTSF) 2 C10 4 conductor. As for (TMTSF) 2 PF 6 
conductor, the question of possible triplet supercon- 
ducting pairing [6,7,10,18-20] is still not completely re- 
solved. Triplet superconductivity is a rather unusual phe- 
nomenon. In our opinion, it has been firmly established 
in the heavy fermion superconductor UPt3 [21] and most 
likely exists in the Sr 2 RuC>4 [22,23] and ferromagnetic 
superconductors [24]. 

In this intriguing situation, it is important that a novel, 
strong candidate for triplet superconductivity - the Q1D 
layered superconductor Lio.gMo60i7 - has been very re- 
cently suggested [25]. As shown in experiments [25], the 
upper critical magnetic field, parallel to a conducting 
axis of this superconductor, is five times larger than the 
so-called Clogston-Chandrasekhar paramagnetic limit for 
singlet superconductivity [26]. A distinctive feature of 
the measurements in Ref. [25] is that the paramagnetic 
limit for the superconductor Lio.gMo 6 Oi7, H p = 3.1 T, 
has been extracted from direct measurements of the 
Pauli susceptibility and specific heat jump (see also 
Refs. [27,28]). Therefore, it does not depend on the de- 
tails of a theoretical model. As also shown in Ref. 
[25], the above-mentioned superconductor is in the clean 
limit and, thus, the spin-orbital scattering cannot be re- 
sponsible for the extremely large experimental value of 



the upper critical magnetic field along conducting axis, 
H* 2 — 15 T. These facts would posit strong arguments in 
favor of the existence of a triplet superconducting phase 
in the Lio.gMogOiy, making it insensitive to decoupling 
Pauli paramagnetic effects. However, as mentioned in 
Ref. [25], orbital destructive effects are minimized when 
the magnetic field is parallel to a conducting axis of a 
Q1D superconductor. Therefore, it is stressed in Ref. [25] 
that the experimentally observed destruction of super- 
conductivity at H > H*2 — 15 T can only be ascribed to 
Pauli paramagnetic effects, which is not in favor of the 
above mentioned triplet scenario of superconductivity. 

The goal of our Letter is to show theoretically that the 
orbital effects can destroy superconductivity even for a 
magnetic field, applied parallel to a conducting axis of 
a Q1D layered superconductor, provided that the inter- 
plane distance is less than the corresponding coherence 
length. By extracting electronic band and superconduct- 
ing phase parameters from measurements in Ref. [25] , we 
show that the above condition holds for the Q1D layered 
superconductor L10.9M06O17. We consider our results 
as the second (after Ref. [25]) major step in establishing 
triplet superconductivity in the above discussed super- 
conductor. 

Let us consider a tight-binding model for the electron 
spectrum of a Q1D layered conductor, 

e(p) = -2t x cos(p x a x ) - 2t y cos(p y a y ) - 2t z cos(p z a z ), 

(1) 

where t x 3> t y 3> t z are overlap integrals of the electron 
wave functions along x, y, and z crystallographic axes, 
respectively. In a magnetic field, 

H=(ff,0,0), A = (0, 0, Hy), (2) 

parallel to the conducting chains of the Q1D layered con- 
ductor (1), it is convenient to write electron wave func- 
tion with definite energy and momentum component p x 
in the following way: 

tptpa:(- x >y> z ) = exp(±ip x x)exp[±ip y z (p x )y] <f>± Px (y,z). 

(3) 
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Note that, in Eq.(3), +(-) corresponds to the left (right) 
sheet of the Q1D Fermi surface (FS), and the functions 
PyiPx) are defined by the following equations: 

vf(Px TPf) T 2t y cos\py(p x )a y ] = 0, (4) 

where vf and p F are the Fermi velocity and Fermi mo- 
mentum, respectively. In this case, we can linearize Eq. 
(1) near two sheets of a Q1D FS in the following way: 

<5 e± (p) = ±2t y a v [p y -Py (p x )]sin[p± (p x )a y ] 

-2t z cos(p z a z ), (5) 

where electron energy, 5e = e — e Fl is counted from the 
Fermi energy e F . 

In a magnetic field, we use the Peierls substitution 
method for Eq.(5), p y - p y (p x ) -» -id/dy, p z a z -> 
p z a z — uo z y/v F , where lo z — eHv F a z / c, e is the elec- 
tron charge, and c is the velocity of light. As a result, 
we obtain the following Schrodinger-like equation for the 
electron wave functions in the mixed (y,p z ) representa- 
tions: 

^Tiv y \p y (p x )]-^ - 2t z cos(p z a z - ^-yj-2fi B sH 

x 0i P Jy,Pz) = 8e 4>t tVx {y,p z ),{6) 

with s being the projection of an electron spin on 
x axis; /ib is the Bohr magneton, Vy[p^{Px)] = 
2t y a y sm\p^{p x )a y }. Note that Eq.(6) can be solved ex- 
actly: 



<Pt, Px (y,Pz) = exp 



x exp 



ztidey 
v y \p y (p x ) 



exp 



±2ifiBsHy 

v y[pf(Px)} 



2t 2 



f 

JO 



cos\ p z a z u)duf. (7) 



It is important that the finite temperatures Green func- 
tions for the wave functions (7), (3) can be determined by 
the standard equation [29]: 

/ + OC 
d(5e)[i>t, Px (xi,yi,p z )]* 
-oo 

x*Pt,pA x ,y,Pz)/(iu n - Se), (8) 

where oo n is the so-called Matsubara frequency. 

In this Letter, we consider the simplest triplet scenario 
of superconductivity in the Li .gMo 6 Oi7, where super- 
conducting pairing is not sensitive to Pauli paramagnetic 
effects [7]: 



MPx,y) = I sgn(p x ) A(y), 



(9) 



where I is a unit matrix in spin-space, sgn(p x ) changes 
the sign of a triplet superconducting order parameter on 
two slightly corrugated sheets of the Q1D FS, A(y) takes 
into account the orbital destructive effects against super- 
conductivity in a magnetic field. It is important that 



the triplet order parameter (9) corresponds to a fully 
gapped Q1D FS (4), which is in qualitative agreement 
with the experimentally observed large specific heat jump 
in Li .9Mo 6 Oi7 superconductor (see Ref.[27]). To derive 
the gap equation for superconducting order parameter 
A(y), we use Gor'kov's equations for unconventional su- 
perconductivity [30-32]. As a result of the calculations, 
we obtain: 



ate) = ■'(/,- 



2irTdyi 



vi\>- 



A(yi) 



, 8t z v F 
x Joi -, — r sin 



u z {y-yi) 

2v F 



sin 



u z (y + yi) 

2v F 



Py 



where < ... > p stands for the averaging procedure over 
momentum p y , g is the electron coupling constant, and 
il is the cutoff energy. 

Note that Eq.(10) is very general. For instance, at high 
magnetic fields and/or low temperatures, it describes the 
exotic reentrant superconducting phase, introduced for a 
different direction of the magnetic field in Ref. [7] . Anal- 
ysis of Eq.(10) shows that we can disregard the reentrant 
superconductivity effects at high enough temperatures, 



T > T*{H) ~ 



2it 2 vf 

and low enough magnetic fields, 

8t z v F 



oj z {H) < 



,,o ' 



(11) 



(12) 



where = 2t y a y . It is possible to show that Eq.(10) 
can be rewritten under the conditions (11), (12) in the 
following way: 



A(y) = g 
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J\y- 
°l v v(Py 



2itTd yi 



y(Py) 



-Vi) 
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— — v y (p y )sma[ v ^ p ^ 
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2v F 



Myi) , (13) 



Py 



It is important that Eq.(13) is still rather general. In 
fact, it describes both the so-called Lawrence-Doniach 
(LD) model [33,34] and anisotropic 3D superconductiv- 
ity. As mentioned in Ref. [25], the LD model condi- 
tion, £ z -C a z /\/2, is not obeyed in the superconduc- 
tor Li .gMo 6 Oi7 and, as we show below, it is possible to 
simplify Eq.(13) to describe anisotropic 3D superconduc- 
tivity: 



= g 



2irTdyi 



\y-yi\>- 



X Jn 



2t z uj z (y 2 - y\) 



Vy(Py)vF 



A(j/l)> • (14) 



Py 
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In Eq.(14), it is convenient to perform the following 
transformation of the variable y\\ y\ — y — zv y (p y )/vF- 
As a result, Eq.(14) can be rewritten as 



My) 



2nTdz 



- ^sinh^] 
'2t z cu 



X Jn 



2y 



V 

v v(Py) 



v v(Py). 



vf 



(15) 



It is possible to show that in the vicinity of supercon- 
ducting transition temperature, (T c — T) <C T c , Eq.(15) 
leads to the Ginzburg-Landau (GL) formula for the up- 
per critical held, 



HUT) 



4ir 2 chT 2 



7((3)et y t z a y a z 



T c -T 



(16) 



where C(x) is the Riemann zeta function. We note that 
the GL slopes for the magnetic field applied perpendic- 
ular to the conducting axis of a Q1D layered supercon- 
ductor were derived in Rcfs. [35,36]. Using the results of 
Ref . [36] , we can write 



HUT) = 



±V2n 2 cT 2 
7((3)ev F t z a z 
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yU,y 



(17) 



(18) 



where the GL coherence lengths along x, y, and z crys- 
tallographical axes are 



7C(3)» 



16(ttT c 
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7C(3)tgag 
8(nT c ) 2 
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8(itT c ) 2 ' 



(19) 



From Eqs.(16)-(19) and experimental data [25], 
fl£(0) ~ 22 T, H y c2 (0) ~ 4 T, ff^(O) ~ 1 T, it is 
possible to estimate the parameters of the Q1D layered 
electron spectrum (1),(4),(5) and the coherence lengths 
(19). They are summarized in Table 1. 
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By using data from Table 1 and Eqs.(ll),(12), we 
can estimate the region of temperatures and magnetic 
fields where Eq.(13) is valid. As a result, we obtain 
T > T* ~ 0.06 K and H < 300 T - conditions, 
which arc well satisfied in experiment [25]. As already 
mentioned in Ref. [25] and as seen from Table 1, in the 
Li .gMo6Oi7 superconductor £ z ~ 2QA > a z /^/2 = 6.7 A. 
Thus, the LD model is not applicable and we can use 




FIG. 1: Solution of Eq.(15) at T = 0.1 K is an oscilatory 
function of the coordinate y. The solution is normalized by 
the following condition: /_ °° A 2 (y)dy — 1. 




T(K) 



FIG. 2: Temperature dependence of the upper critical mag- 
netic field, Hc 2 (T), numerically calculated from Eq.(15) al low 
enough temperatures, is shown by a solid line. The Ginzburg- 
Landau linear dependence (16), which is valid at T c — T <C T c , 
is shown by broken line. Rectangles represent the experimen- 
tal data [25]. 



Eqs.(14),(15) for anisotropic 3D superconductivity. We 
note that Eqs.(14),(15) are qualitatively different from 
the gap equations for a 3D isotropic case [37,38], since 
the former take into account Q1D topology of the FS 
(4), (5). In particular, a typical solution of Eq.(15) at low 
temperatures, T <C T c ~ 2.2 K, changes its sign with 
changing coordinate y, in contrast to 3D isotropic case, 
as shown in Fig.l. We solve Eq.(15) numerically in the 
range of temperatures, T c > T » T*, and compare the 
obtained temperature dependence of the upper critical 
magnetic field, H* 2 {T), with the experimental data [25] 
in Fig. 2. As seen from Eq.(2), the calculated dependence 
of H^ 2 (T) is in good qualitative and quantitative agree- 
ment with the experimental results. 

To summarize, we have explained theoretically the ob- 
served destruction of superconductivity in the Q1D lay- 
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ered superconductor Lio.gMo60i7 in a magnetic field par- 
allel to its conducting axis [25], in the framework of 
triplet superconductivity scenario. We have also sug- 
gested the most probable triplet order parameter [see 
Eq.(9)]. It corresponds to the absence of Pauli param- 
agnetic destructive effects against superconductivity and 
is qualitatively consistent with a large value of the ex- 
perimentally observed specific heat jump at the super- 
conducting transition at H — [27]. [Note that from the 
microscopic point of view, triplet phase can be stabilized 
in a Q1D superconductor as a result of repulsive inter- 
chain electron-electron interactions [39-41].] 

Below, we would like to discuss the applicability of the 
Fermi liquid (FL) approach we have used in the Letter 
to describe the superconducting phase transition in the 
Lio.gMogOiy. First of all, we note that at high enough 
temperatures, the Luttinger liquid effects are observed 
[42-44] in the Lio.gMogOiv conductor. This naturally re- 
flects the Q1D nature of its electron spectrum (4), (5) 
and is not crucial for our analysis, since we consider the 
low temperature region, T < T c = 2.2 K. In this con- 
text, it is important that from the theoretical point of 
view, the FL picture is restored at temperatures lower 



than t z , t y 



10 — 45 K. Another point of concern 



is the experimentally observed increase of resistivity at 
T < T min ~ 15 — 30 K. So far, its nature has not been 
clearly understood. At present, there exist two most 
popular competing points of view on this resistivity in- 
crease phenomenon: localization effects [45] , and the pos- 
sible partial charge-density-wave instability or the cor- 
responding fluctuations [46,47]. Nevertheless, we stress 
that there are two experimental features [25], which are 
important for the validity of our analysis. The first one 
is based on the fact that the noted increase in resistivity 
is of the order of Sp/p ~ 0.25 [25] and, thus, small. The 
second feature is that the magnetoresistance at T ~ 4 K, 
as shown in Ref. [25], demonstrates quadratic behavior, 
5p(H)/p(0) ~ H 2 , which is a direct test of the FL theory. 
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